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About the connection between the Ct power spectrum of the Cosmic 
Microwave Background and the T m Fourier spectrum of rings on 
the sky 



m 
O 
O 

(N 

< 

OO 

m 
> 

in 

(N 

o 
m 
o 



6 



X 



R. Ansari 1 , S. Bargot 1 , A. Bourrachot 1 , F. Couchot 1 , J. Haissinski 1 , S. Henrot-Versille 1 , 
G. Le Meur 1 , O. Perdereau 1 , M. Piat 2 , S. Plaszczynski 1 and F. Touze 1 

1 Laboratoire de VAccelerateur Lindane, IN2P3-CNRS and Universite de Paris-Slid, BP 34, 91898 Orsay Cedex, France 

2 Institut d Astrophysique Spatiale, INSU-CNRS and Universite de Paris-Slid, 91405 Orsay Cedex, France 

original form 2003 January 16 

ABSTRACT 

In this article we present and study a scaling law of the mT m CMB Fourier spectrum on 
rings which allows us (i) to combine spectra corresponding to different colatitude angles (e.g. 
several detectors at the focal plane of a telescope), and (ii) to recover the Ce power spectrum 
once the T m coefficients have been measured. This recovery is performed numerically below 
the 1 % level for colatitudes O > 80° degrees. In addition, taking advantage of the smoothness 
of the Ce and of the T m , we provide analytical expressions which allow the recovery of one 
of the spectra at the 1 % level, the other one being known. 
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1 FOURIER ANALYSIS OF CIRCLES ON THE SKY 
VERSUS SPHERICAL HARMONICS EXPANSION 

Cosmological Microwave Background (CMB) exploration has re- 
cently made great progress thanks to b allo on born e xperiments 
(BOOMERANG 2000, MAXIMaI EoOoI and lARCHEOPSiEool 
ground-based interferometers IcBllboOa iDASjllOOZlVSAl 
MAP 1 whose first results will be available at the begin- 
ning of 2003 and the forthcoming Planck satellite 2 whose launch 
is scheduled for the beginning of 2007 will scan the entire sky with 
resolutions of 20 and 5 minutes of arc respectively. These CMB ob- 
servation programs yield to a large amount of data whose reduction 
is usually performed through a map-making process and then by 
expanding the temperature inhomogeneities on the spherical har- 
monics basis: 



(ARCHEOPS, Map, Planck) perform or will perform circular scans 
on the sky. Carrying out a one-dimensional analysis of the CMB in- 
homogeneities on rings provides a valuable alternative to character- 
ize its statistical properties (Del abrouille et alll998l) . A ring-based 
analysis looks promising e.g. for the Planck experiment where re- 
peated (~ 60 times) scans of large circles with a colatitude angle 
~ 85° are being planned. This approach differs in several ways 
from the one based on spherical harmonics. In particular it does not 
require the construction of sky maps and some systematic effects 
could be easier to treat in the time domain rather than in the two- 
dimensional (0, tp) space (1// noise for instance) since the map- 
making procedure involves a complex projection onto this space. 
For a circle of colatitude 0, one writes 
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The outcome of the measurements is given in the form of the an- 
gular power spectrum Ce =< \ai m \ >. The set of Ce coeffi- 
cients completely characterizes the CMB anisotropie s in the case 
of un correlated Gaussian inhomogeneities iHu & Dodeisonll200l 
iBond & EfstathiovJl987h . 

Several of the current or planned CMB experiments 



1 Map home page: 

|http : / /map . gsf c . nasa.gov/| 

2 Planck home page: 

|http ://astro.estec.esa. nl/SA-general/Pro jects/Planck/| 



and the r m Fourier spectrum is defined by 



(2) 



(3) 



These T m coefficients are thus specific to a particular colatitude 
angle 0. We propose below a simple way of combining sets of 
such coefficients corresponding to different values (i.e. different 
detectors). 

Fig. Q shows an example of the Ce power spectrum for £ < 
1500, together with two Fourier spectra 3 which describe the same 



3 Note that we have chosen the following normalizations: the Ce coeffi- 
cients have been multiplied by £(2£ + 1)/4tt and the Y m by Ira. 
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Figure 1. 'Typical' power spectra. The insert shows a l(2l+l)Ct / in spec- 
trum up to I = 1500. The main graphs are two Fourier spectra (2mT m ) 
exactly calculated using Eq.|4]: one for = 90° (darker curve), and the 
other for © = 40° (lighter curve). The triangles represent a subsample of 
the 2r«r m (© = 40°) coefficients after having rescaled their abscissa by a 
factor 1/ sin 40° = 1.556. 



Figure 2. Absolute differences (in [iK 2 ) between 2mr m (6) spectra 
scaled to = 90° and the interpolated 2mr m (© = 90°) spectrum. All 
spectra are based on the Ci spectrum of FigQ We have worked out these 
differences for = 60° (smallest amplitude curve), 40°, 30°, 20° and 
10° (largest amplitude curve). The insert displays the low m part, showing 
that the difference has a meaningful value only above m = 2/ sin 0. 



sky for two quite distinct cases, one for = 90° and one for = 

40°. 

Note that for this Fig.Qand throughout the article the Co and 
Ci coefficients have been set equal to 0. 

The rela tion that gives the r m (0) from the Cg was obtained 
bv lDelabrouille et aljjl993l : 

oo 

r m (e)= Yl c e B 2 e rL(cose), (4) 

t=\m\ 

where the set of Be coefficients characterizes the beam function and 
the V\ m are the normalized associated Legendre's functions. This 
relation assumes that the ae m introduced in Eq.Qare uncorrelated 
Gaussian random variables and that the scan is performed with a 
symmetric beam. 

In this article, we present the scaling law and the inverse trans- 
formation that consists in the calculation of the Ci from the Y m . In 
section 2, we demonstrate that this simple scaling law, displayed 
by the mY m spectrum for different colatitude angles, is accurate. 
Section 3 is dedicated to the description of two different methods 
proposed to invert Eq.|4]in the case O = 90°. While a simple ma- 
trix inversion leads to the result, we also present an approximate 
analytic method. In section 4 these two methods are extended to 
the general case where < 90° . 



2 SCALING OF THE mT m (©) SPECTRUM 

Our study was triggered by one of us noticing that the product 
mV m (Q) is only function of the reduced variable /i = ra/sinO, 
i.e. this product is independent (to a very good approximation) of 
the colatitude angle O. 

This scaling is illustrated in Fig. Q where a 2mV m spectrum 
computed for a colatitude angle of = 40° is scaled to match 



the corresponding = 90° one. To quantify the precision of 
this approximate scaling law, we have computed the differences 
between the scaled 2mF m (0) and the interpolated 2mF m (0 = 
90°) spectrum (at m/ sin0). Examples are shown in Fig. [2] for 
five values ranging between 60 and 10° . The absolute values of 
these differences are lower than 2 fj,K 2 over the whole m range for 
the particular spectrum given in figureQ They are only defined for 
m values greater than 2/ sin 0, as shown in the insert. Oscillations 
are observed in the difference. They present the same period but 
their amplitudes increase as the colatitude angle decreases. 

Different 2mr m (0) sets obtained from several detectors over 
a small range of colatitude angles (a few degrees) may be com- 
bined using this scaling law, with a precision better than 0.01% . 
Several experiments, spanning a wider range of colatitude angles, 
may also be combined likewise, however with a slightly worse pre- 
cision. 

In the following, we explain this scaling law using a geomet- 
rical and a mathematical argument. 

2.1 Geometric interpretation 

The power spectrum r m (0) is the Fourier transform of the 
signal autocorrelation function A(5(/>, 0), where 8cj> is the phase 
difference between two points of the scanned ring. Two such points 
have an angular separation Sip on the unit sphere, where: 

Sip = 2arcsin(sin0 sin ^) . (5) 

This relation between 8<f> and 8ip allows one to express the 
scaling law, since the signal autocorrelation function, expressed as 
a function of Sif) is equal to the autocorrelation function on a large 
circle scan: 

A(5ip,ir/2) = A(6 <!>,&). (6) 
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For small S</>, this relation becomes linear: 

Sip = sin 6 6<t> . (7) 

So that, in this linear regime, the autocorrelation function sat- 
isfies: 

A{6cf>, G) = A(sin 6 5(f), tv/2) . (8) 

Since the ring length L on the unit sphere is 2ir sin O, the m th 
harmonic of the Fourier expansion corresponds to structures on the 
sky of angular size 

sin O 2ir 



A = 2tt- 



(9) 



In the continuum approximation, taking the Fourier transform 
of both sides of Eq.[5]leads to: 



r m (e) = 



ml sin 



sin O 

which, using Eq. |5|leads to the scaling law: 
mr m (6) =^r M (7r/2). 



(10) 



(11) 



While we are mainly concerned here with circular scanning, 
the same reasoning can be made for any kind of trajectory on the 
sky as long as it stays 'close' to a large circle on angular scales 
of order A, and the same scaling law applies to the power density 
spectrum expressed as a function of 1/A. 



2.2 Analytic interpretation 

To investigate this scaling mathematically, we start from Eq. |4| 
which gives the exact relations that connect the r m (0) to the Ci. 
Since the Be are - supposedly - well known quantities for each 
experimental set up, we will no longer mention them explicitly and 
we will deal with the coefficients Ce = GgBf. 

We calculate the "Pf m (cos O) factors using approxim ate ex- 
pressio ns of the Legendre's associated functions given by iRobir] 
( 1957) (see appendix A for some details) which, once normalized, 
read 

• for £ < raj sin : 



, 1 £+\ (£cosQ + M\ e+ i 
(cos 6) ~ — \ 



m cos Q — M\ m I £ + k — m 



(I — to) sin O 



n 



£ + k 



(12a) 



where M = \J m? — £ 2 sin 2 0, 
• and for £ > raj sin : 



Pfm (cos 0) 



2tt 



(-l) m / 2(22 + 1) 



N 



(126) 



where N = yl 2 sin 2 — m 2 . The expression of the angle u is 
given in appendix A. These approximations are illustrated by Fig.|3] 
For £ < to/ sin 0, the numerical value of Pf m (cos 0) is neg- 
ligible, while for £ > to/ sin Eq. 6b implies 



P £m (cos0) 



2£ + 1 



4tt 2 (£2 gin 2 - to 2 ) 1 /: 



-[l + cos(2u)]. (13) 



Since the CMB angular power spectrum varies slowly as a 
function of I, we may replace the sum over £ in Eq.|4|by an integral. 
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Figure 3. Comparison between the exact value of Vt 80(>( COS (40°)) as 
a function of i, dotted line, and the one obtained with the approximate 
expressions of Eqs. (12a) and (12b), solid line. The arrow indicates the 
£ = 800/ sin 40° abscissa. 



We thus obtain 



C(£)[2£ + 1][1 + cos(2lo)} 



(£ 2 sin 2 - m- 



1)1/2 



(14) 



where £ ma x is an £ value beyond which the power spectrum van- 
ishes, and C(£) is a function of £ £ ]0, £ ma x] that smoothly inter- 
polates the Ce coefficients (a simple way of proceeding is given in 
Appendix B). 

The oscillation frequency v of the cosine term (as a function 
of £) in the integrand in the right side of Eq. ll4l is of order 0/7T 
(thus v ~ 1/2 when = 7r/2). Such a frequency is high enough 
for this cosine term to contribute to a very small amount to the 
integral. This will be checked numerically in section |3~T1 below. 
Thus we may write: 



mr m (0) 



4tt 2 



C(£) 



21+1 



[{£/nY - l]i/: 



■dl 



(15) 



This equation demonstrates - within the approximations that have 
been made - that the product mr m (0) depends only on the vari- 
able fi = m/sin 0. 

Since the variable [i is not constrained to be an integer, one 
has to introduce a smooth function, r(m, 0) where m is now a 
real, that interpolates the T m (0) discrete spectrum. This can be 
done in the same way as the one indicated for the Ci spectrum (cf. 
Appendix B). 

In terms of this T(m, 0) function, the scaling law is expressed 
by the relation: 



tv™' c«"i sin6) TV™' sine 



(16) 



sin 0' sin 0' 

This equation follows from the equality mr(m, 0) = 

m'r(m',0') which holds true provided that m/sin0 = 
to'/ sin©'. 

Assuming that the Fourier spectrum has been obtained for 

a particular value of the colatitude angle, Eq. 1161 allows one 

to calculate r(m',0') for ml — msin0'/sin0, m = 
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Figure 4. Comparison between the 2mT m coefficients computed with the 
associated Legendre's polynomials (dashed line) and the 2raT(m) function 
calculated using Eqs. ll7llBT1 and lB2l with a = 0.5 (solid line). The relative 
difference between the results of the two calculations is shown by the lower 
curve (in %, right scale). 




Figure 5. Comparison between the 'Typical' Ci coefficients for © = 90° 
(solid line), used to calculate the T m Fourier spectrum (using the P(0) 
matrix) and the C(t) function obtained by inserting this Fourier spectrum 
in Eg. 1191 (triangles; only some points are shown). We have set a = 1 in 
Eq, IB2l The relative difference (in %) is shown by the lower curve (right 
scale). Insert: zoom on the low I region. 



1,2 •■•m max — ^maxSinO. Then, by interpolation, one gets 
r(m', O') for all integer values of m! ranging from sin 0'/ sin Q 
up to £ m ax sin B'. Eq. ll6l can thus be used to compare and combine 
Fourier spectra that correspond to different O values. 



3 RECOVERING THE C t COEFFICIENTS FROM THE 
r m 0/2) FOURIER SPECTRUM 

3.1 Checking and solving the integral equation that relates 

C(£) tor(m,7r/2) 

Since O is assumed to be equal to 7r/2 in this section, the variable 
[i can be identified with m. 

In order to facilitate the numerical calculation of the right side 
of Eq.[l5] we introduce a new variable of integration x defined by 
£ — m cosh x. Then Eq. ll5l can be rewritten 



4tt 2 



COsh (&max / m) 

(2m cosh x + l)C(m cosh x) dx . 



(17) 

The transformation defined by Eq. ll7l is linear: thus one may 
insert in the integrand an interpolating function of the d spectrum 
as defined by Eq. lBll The output of Eq. ll7l applied to the angular 
power spectrum of Fig. Q is shown in Fig. [4] One can see in this 
Fig.|4|that for such a spectrum the approximations that were made 
in section[2]ensure an accuracy better than 1% - except at the lower 
end of the spectrum where the relative error drops below 2 % for 
m = 14. 

Eo. I17l can be solved for C(£) by noticing that this integral 
equation is similar to Schlomilch's equation which reads 



F(m) 



2 

i" JO 
where m is a real 



/2 



$(msinx) dx 



(18) 



The way to solve the latter equation can be found e.g. in 
lKraznov|i 19771) . We proceed in a similar way for Eq. ll7l (the details 
are given in Appendix C) and we obtain 



C{£) = -8tt 7 



cosh- 1 (£ mQX /«) 

ofj _ | ] r'(^coshx) dx, (19) 

where V' is the derivative of T(m, n/2) with respect to m. Again 
the transformation implied by Eq. ll9l is a linear one, allowing the 
use of interpolating functions as defined in Appendix B. Fig.[5]il- 
lustrates the use of this integral equation to calculate the Ci coeffi- 
cients starting with the set of r m (7r/2)'s. 



3.2 Numerical inversion 

In the 6 = n/2 case, the connection between the set of C/s and 
the correspondi ng T m 's is si mple since Eq.|4]can be written using 
matrices I M. Piat et al. 20Q3): 



r = p(o) xc, 

with: 

P(0)„ = [Pji(0)f 



(20) 



(21) 



where Vji are the normalized associated Legendre's functions. 
P(0) is (upper) triangular. 

In addition, since the associated Legendre polynomials are de- 
fined as: 



Pim(0) = 



(-1) 





p (2<+2m)! 



[a], 



n , r? — 77 if I — m — 2p 

2'p!(p+m)! r -" (22) 

if£-m = 2p + l [b], 



all of the P(0).. diagonal elements are different from zero - thus 
this matrix is invertible. 
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Figure 6. Middle curve (solid line) : Fourier spectrum obtained using 
Eq. ll7l when only C300 7^ 0. We have used a = 0.5 in Eq. lB2l Upper and 
lower set of points: the T m coefficients computed with the P(0) matrix. 
Since we assume here that = it/2, all T m coefficients whose indexes m 
are odd vanish. 



The inverse of P(0) is also upper triangular and keeps the 
peculiar structure of the original matrix: in both P(0) and P(0) _1 
only the I — m = 2p terms differ from zero. 



3.3 Comparison between the analytic and the numerical 
transformations 

One way of comparing the two methods of calculating the Fourier 
spectrum is to look at what happens when a single Ce coefficient 
is different from zero. This is done in Fig. [6| for the case where 
C300 = 1- Note that since we assume here that O = 7r/2, Eq. 22 
b implies that all T m coefficients with an odd index vanish (for a 
single non vanishing Gt coefficient with an odd I value, all T m co- 
efficients with an even index would vanish). One notices that the 
r(m) function runs at mid-height of the non-vanishing T m coeffi- 
cients. 

Conversely, one may look at the C(£) function that corre- 
sponds to the case where a single F m Fourier coefficient is differ- 
ent from zero as shown in Fig.Q(here we used F^oo 7^ 0). The fact 
that the C(l) graph is negative in some domain of I values shows 
that no distribution of temperature inhomogeneities which satisfies 
the validity conditions of Eq. [4] (isotropy and Gaussian ai m ) can 
correspond to a Fourier spectrum with a single non- vanishing coef- 
ficient. 

Taken together, Figs. [6] and [7] show where we should expect 
a strong signal in one spectrum in the case the other spectrum 
presents a high power in some particular bins. 



Figure 7. Solid line: the C( spectrum obtained with Eq. 1191 when only 
F300 7^ (we have set a = 1 in Eg. IB 21 . In dots: the Ci coefficients 
calculated with the P(0) _1 matrix . 



4 WORKING WITH SMALLER RINGS ON THE SKY 

(0 < tt/2) 

4.1 General features of the Fourier spectrum 

In the preceding section we assumed that the scanned rings are the 
largest ones on the sphere (O = 7r/2). In this case the fact that the 
P(0) matrix is invertible establishes that the Fourier spectrum of 
such rings contains all the physical information carried by the Ci 
coefficients. 

Scanning smaller circles on the sky implies a higher funda- 
mental frequency in Fourier angular space and thus a less dense 
sampling of this Fourier space. 

In fact the loss of information is then twofold: 

• The G(fi) = mT{m, O) function is no longer measured for 
jti = 1: the lowest value of /j, which can be reached with the data is 
now jj, — 1 /sin O. 

• Secondly, G(fi) is no more measured for [i values that differ 
by one but for [i values that differ by 1/sinO. As a very simple 
example: if the scan is performed for G = it/6, then one mea- 
sures G(fj.) only for fi — 2n with n G ]0, £mam/2\- Because of the 
smoothness of the angular spectra, this sparse sampling of the func- 
tion G(n) is not necessarily a drawback as long as the accuracy of 
the measurements compensates for it. 



4.2 Analytic calculation of the Ce spectrum for L > 1/ sin 

As far as the analytic calculation of the Ct spectrum is concerned, 
it can be performed with the same formalism as above (cf. subsec- 
tion l3. li . One should merely replace the derivative of T(m, 7r/2) 
that appears in the right side of Eg. ll9l bv the derivative (with re- 
spect to m) of 

imax Sine 

f(m)=sine r</(msin0-i). (23) 
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f (m) is just the rescaled version (cf. Ea.rnfl of r(m, 8) denned 
by Ea, IB3l (this rescaling translates the F(m, G) Fourier spectrum 
into the one corresponding to G = 7r/2). Furthermore the width of 
the interpolating function f(x) of Appendix B (see Eq. lB2t should 
be increased by a factor 1/sin B. 



4.3 Numerical calculation of the Ct spectrum for £ > 1/ sin G 

It follows from subsection 14,11 above that the r m (G) coeffi- 
cients differ significantly from zero in the range 1 ^ m ^ 
(■max sin B. Then using the F(m, 6) function that interpolates 
these coefficients and Eg. H6l one can calculate the following set 

nf / — 



x — t-min + 1 values 
sin B r(m' sin 6, 6), 



(24) 



with m! = t 



+ 1,- 



tmax where £ min is the first integer 
larger than 1/sinO. These F m i coefficients are the ones of the 
Fourier spectrum for O = 7r/2. Once obtained, the Ct spectrum is 
simply given by 



c = p(o) _1 r 



(25) 



fori 1 ^ Imin - The P(0) matrix and its inverse have been discussed 
in section [3T2l The first £ m in ~ 1 rows and columns of P(0) -1 
should be omitted in Eq. l25l since the lowest value of the m' index 

is £min • 

Fig. [H shows a numerical example: we use the 'typical' Ct 
spectrum of Fig.Q to produce a set of T m values in the = 40° 
case (Eq.|4}- Then we apply the method described above and com- 
pare the input spectrum to the obtained one. In this example we 
used a simple linear interpolation of the T m spectrum. The agree- 
ment is excellent and better than the one obtained with the analytic 
method (cf. Fig.|5J as the latter involves some approximations (cf. 
section 2) in addition to the ones stemming from the scaling and 
the interpolation procedure. 

The excellent agreement of Fig.[S]breaks down for low values 
of I. Nevertheless, for £ > 3 in the O = 40° case, one gets an 
agreement better than 10% (far above the cosmic variance). For the 
case G = 80° our simple scaling method can be used up to an 
accuracy better than 1% for any i values. 



5 CONCLUSION 

We have shown how data taken on circles with different colatitude 
angle 6 can be combined using a scaling law that is satisfied by the 
mT m (Q) coefficients at the 0.1 % level in a wide range of m and 
G values. 

Then we have derived this scaling property from both geomet- 
rical considerations and linear expressions of the T m coefficients 
in terms of the Ct ones by introducing analytic approximations 
of the normalized Legendre's associated polynomials Vt m (cos G) 
that enter these relations. 

Integral equations were obtained that relate to a good approx- 
imation interpolating functions of the two sets of coefficients (F m 
and Ct). These analytic relations give a simple picture of the con- 
nection between the two types of spectra and are easy to use. 

Finally we have investigated ways of calculating the Ct coef- 
ficients when the Y m Fourier spectrum is known. We have shown 
how the inverse of the P| m (0) matrix can be used to perform this 
calculation not only for G = 7r/2 but also in the general case where 
G < 7r/2. This was achieved by taking advantage of the scaling of 



1(21+1)C /4ji (p.K 2 ) 



Relative 
difference ' 
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Figure 8. The input Ct spectrum (solid curve) and the one reconstructed by 
the numerical method in the = 40° case (only some points are shown). 
The relative difference between the two spectra is shown by the lower curve 
(in %, right scale). Insert : zoom on the low £ region. 



the mF m spectrum on the one hand and of its smoothness on the 
other. 

This set of results provides a basis for further investigation of 
the connection between the measured Ct and F m spectra altered 
by noise and errors. 



APPENDIX A: APPROXIMATE EXPRESSIONS OF THE 
NORMALIZED LEGENDRE'S ASSOCIATED 
POLYNOMIALS 

We start with asymp t otic ex pressions of the Legendre's functions 
obtained by iRobinl 1 19571) in the limit of large I, raj I being 
kept constant. These asymptotic expressions depend on the relative 
value of m and £ sin G. 
• For £ < m/ sin G, 



Pf m (cosG) 
(~l) m £\ 



t^ 1 

U cos e+M) m cos Q — Af ) 



27r(7-m)! 



-m) m M2 . 



' e 



(Al) 



where M = \J m 2 — P sin 2 G , 
• while for £ > mj sin G, 



(cos G) 



(-i)"V- 

7T 



£!(l-m) 



1 l + r, 

~4(l+ m ) 2 



(£-m)!/ + S JV5 

where 

N = \/£ 2 sin 2 G-m 2 , 

u = (£+-)a-mf3-j, 

a — arg [£ cos Q + iN) , 







arg(m cos G + iN) . 



(A2) 

(A3) 

(A4) 

(A5) 
(A6) 
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To 'normalize' these polynomials and obtain the Vim ones, 
they must be multiplied by 



21 + 1 {£ - m)\ 
4-7T (£ + m)\' 



(A7) 



Then the last step consists in using Stirling's formula 
(n! ~ \^27vn n+ 2 e~ n ) to replace the factorials by analytic 
functions. A few simplifications can then be made that lead to the 
approximate expressions used in section 2. 



Once the integration order is reversed in the right side of this equa- 
tion on obtains: 

poo poo 

/ r'(wcosh^) dip = / h (it cosh £) sinh £d£ 
Jo Jo 

r * cosh x dx 



o \J sinh 2 £ — sinh 2 



2u 



h{u) = / r'(£coshV0 dip. 



(C4) 



The integral over x is simply n/2. Furthermore h(oo) — in our 
case, so that 



(C5) 



APPENDIX B: INTERPOLATING FUNCTIONS OF THE 
DISCRETE POWER SPECTRA 

Since the calculation of the Ci coefficients involves integr als over 
spher ical Bessel ji functions (see, e.g. ISeliak U. & Zaldarriagal 
1996), one may try and use an expression of these functions that 
extends them to non-integer values of j. But here we will adopt a 
much simpler procedure and write 



C(£) = £ C t f(£- 



(Bl) 



where I is now a real whose value ranges between 2 (recall that we 
ignore the dipole term) and £ ma x, and f(x) is a positive, infinitely 
differentiable function (/ 6 C°°), which differs significantly from 
in an range which is of order unity, and whose integral over x 
is unity. In practice we used 



1 — x" 

^7T(7 2<H 



(B2) 



with a ~ 1. 

Similarly, we define an interpolating function for the T m (Q) 
coefficients in the following way: 



sine 

F(m,@)= ]T r 4 /(m-i), 



(B3) 



REFERENCES 

ARCHEOPS: Benoit, A. et al. |astro-ph/0210305| submitted 

Bond, J. R. & Efstathiou, G. 1987, MNRAS, 226, 655 

CBI: Contaldi, C. R., Bond, J. R., Pogosyan, D., et al. preprint 

|astro-ph/0210303| 
Delabrouille, J., Gorski, K. M.,& Hivon, E., 1998, MNRAS , 298, 

445 

DASI: Halverson, N. W., Leith, E. M., Pryke, C, et al. 2002, ApJ 
568, 38 

MAXIMA: Hanany, S., Ade, P., Balbi, A. et al. 2000, ApJ, 545L, 
L5 

Hu, W. & Dodelson,S. 2002, ARA& A, 40, 171 
Kraznov, M et al., Equations integrales, Editions de Moscou, 
1977. 

BOOMERANG: Mauskopf, P. D.,Ade, P. A. R., de Bernardis, P., 

et al. 2000, ApJ, 536L, L59 
M. Piat, G. Lagache, J.P. Bernard, M. Giard, and J.L. Puget, A& 

A 393,359-368,2002 
Robin, L., Fonctions spheriques de Legendre et fonctions sphero'i- 

dales, Gauthier-Villars (3 vol. 1957-1959) 
Seljak, U. & Zaldarriaga M., 1996, A&A , 469, 437 
VSA: Taylor, A. C, Carreira, P., Cleary, K., et al. preprint 

last ro-ph/ 02 05*3811 



where m is a real and f(x) is chosen as above. 



APPENDIX C: INVERSION OF THE INTEGRAL 
EQUATION RELATING C(£) TO V(m) 

Since C(£) vanishes for £ > £ ma x, the integral equation I17i is of 
the form 

p oo 

r(m) = / h(m cosh x) dx . (CI) 
Jo 

We differentiate both sides of this equation with respect to m, sub- 
stitute for this variable m the product u cosh tp, and integrate both 
sides over ip between the limits and oo. We thus obtain: 

oo p oo r oo 

r^ucoshf/)) dip =1 dx / h! (u cosh tp cosh x) cosh x dtp. (C2) 
Jo Jo 

Then a new integration variable £ is used in the second in- 
tegral of the right side of this equation, defined by cosh^ = 
cosh ip cosh x. Some simple algebra then leads to 

f%, , ,\,, f°? f°°ti (u cosh f) sinh £ cosh x ^ 
/ r'{ucoship)dip =/ dx / K ^' s — d£. (C3) 
Jo Jo Jx \f sinh £ — sinh x 



P l800 (cos(40)) 



1.5 



0.5 







■0.5 



l=800/sin(40) 



Approximate expression 
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